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We study the role of the quench temperature Tj in the phase-ordering kinetics of the Ising model 
with single spin flip in d = 2, 3. Equilibrium interfaces are flat at Tj — 0, whereas at T/ > they 
are curved and rough (above the roughening temperature in d = 3). We show, by means of scaling 
arguments and numerical simulations, that this geometrical difference is important for the phase- 
ordering kinetics as well. In particular, while the growth exponent z = 2 of the size of domains 
L{t) ~ t^^^ is unaffected by Tf, other exponents related to the interface geometry take different 
values at Tf — or Tf > 0. For T/ > a crossover phenomenon is observed from an early stage 
where interfaces are still flat and the system behaves as at Tf = 0, to the asymptotic regime with 
curved interfaces characteristic of Tf > 0. Furthermore, it is shown that the roughening length, 
although sub-dominant with respect to L(t), produces appreciable correction to scaling up to very 
long times in d = 2. 

PACS: 05.70.Ln, 75.40.Gb, 05.40.-a 

I. INTRODUCTION 

When a binary system in suddenly quenched from above the critical temperature Tc to a temperature Tf < Tc, 
phase-ordering occurs with formation and growth of domains. After a certain time tsc dynamical scaling sets 
in, characterized by the typical size of ordered regions growing algebraically in time, L{t) ~ t^/'^ . When domains 
are large, their bulk is in quasi-equilibrium in one of the two broken symmetry phases which are characterized by 
finite correlation length and relaxation time teq{Tf) ~ £,^{Tf), while the motion of the boundaries keeps the 

system globally out of equilibrium. At a given time s, therefore, non-equilibrium effects can be detected by looking 
over distances larger than L{s), because in this case one or more interfaces will be observed. On the other hand a 
local observation performed from time t = s onwards can reveal non-equilibrium features only for time separations 
t — s > s, because on these timescales at least one interface has typically passed across the observation region. In 
the other regime, instead, for space separations r <ti L{s) or time separations t — s <^ s, the equilibrium properties of 
the interior of domains are probed. This character of the dynamics induces an additive structure for pair correlation 
functions between local observables. Using the terminology of spin systems, and considering, for simplicity, the spin- 
spin correlation function G(r, i, s) — {ai(t)aj{s)) — {ai{t)){aj{s)), where ai(t) is the value of the spin on site i at time 
t and r is the distance between sites one has 

G{r,t,s) ^ G,tir,t- s) +Gagir,t,s). (1) 

The stationary term Gst describes equilibrium fluctuations inside domains, and decays to zero for distances r 3> ^{Tf) 
and/or time separations t — s ^ teq{Tf). Gag, which contains the out of equilibrium information, is the correlation 
function of interest in the theory of phase ordering, and obeys the scaling form |^ 

Gagir,t,s) = G{r/Lis),t/s). (2) 

Furthermore, the scaling function for large time separation [1] is of the form 

G{r/Lis),t/s)^{t/s)-'/^h[r/Lis)] (3) 

where A is the Fisher-Huse exponent and, in system with sharp interfaces, like the Ising model, the function h{x) 
obeys the Porod law 



1 — h(x) ~ X 



(4) 
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for a; < 1. In general, both the terms of the spUtting ([T]) display a universal character. For the stationary part, 
this is well known from equilibrium statistical mechanics, where the renormalization group allows the classification of 
different systems into universality classes on the basis of few relevant parameters A similar property is believed 
to hold also for the aging term. Universal indices, such as the exponents A, z, or other appearing in different 
quantities, should depend only on a small set of parameters among which the space dimension d, the number of 
order parameter components and the presence of conservation laws in the dynamics. The theoretical foundations of 
this idea, however, are not as robust as for its equilibrium counterpart This is due to the non-perturbative character 
of the dynamical problem. Actually, while in equilibrium an upper critical dimension du exists above which the 
renormalization group (RG) fixed point is Gaussian, allowing the e-expansion for d < djj, there is not an upper 
critical dimension for the dynamical process following a quench below (6|. Although an approach based fully on 
the RG is not available, complementing RG techniques with a physically motivated ansatz, it has been shown Q that, 
for a system of continuous spins described by a time-dependent Ginzburg-Landau (TDGL) equation there exists an 
attractive strong coupling fixed point at T = governing the large scale properties of quenches to every Tf < Tc- 
This result supports the idea of a universal character of the aging term in Eq. ([T]), allows for a definition of non- 
equilibrium universality classes and shows that universal quantities, such as exponents, are the same in the whole low 
temperature phase. Restricting from now on to scalar systems with short range interactions and without dynamical 
conservation of the order parameter [ij , these quantities should depend only on space dimension 8] . This would agree 
with the physical idea that Tf only determines the size £,{Tf) of the thermal island of reversed spins inside the domains, 
described by the stationary term in Eq. ([T]) , leaving unchanged large-scale long-time properties of the interface motion, 
contained in the aging part. Basically, this indicates a unique mechanism governing the non-equilibrium behavior of 
interfaces. Restricting our attention to the exponent z, the Lifschitz-Cahn- Allen theory confirms this idea, since 
it gives z = 2 for every Tf < T^. At the basis of this result is the so called curvature driven mechanism: the existence 
of a surface tension implies a force per unit of domains boundary area proportional to the mean curvature which, in 
turn, is proportional to the inverse of L{t). For purely relaxational dynamics, this readily gives z = 2, independent 
on Tf and on dimensionality. 

These results are all based on continuous models where the usual tools of differential analysis can be used and the 
curvature is a well defined object. This approach is justified also for lattice models, such as the nearest neighbor 
Ising model, at relatively high Tj, where temperature fluctuations produce soft interfaces, which at a coarse-grained 
level have a continuous character, and can be well described in terms of partial differential equations. When the 
temperature is lowered, however, these interfaces become faceted. This means that, although the growing structure 
has still a bicontinuous interconnected morphology, interfaces are flat up to scales of order L{t). This implies that their 
description in the continuum may be inappropriate. Then, while continuum theories predict temperature to be an 
irrelevant parameter, with Ty-independent exponents and a common kinetic mechanism for all quenches to Tf < Tc, 
lattice models could in principle behave differently, in particular at = 0. This would imply that temperature 
fluctuation do play a significant role in the way interfaces evolve, determining, besides the properties of the stationary 
term in Eq. ([T]), also those of the aging contribution. This issue is not yet clarified; let us mention, for example, that 
while for quenches to T/ = in d = 3 the exponent z = 2 has been observed [10] in numerical simulations of the 
TDGL equation, for the Ising model one measures (T]| an higher value whose origin is not yet clear. 

In this Paper we consider the role of Tf in the phase-ordering kinetics of the nearest neighbor Ising model. For 
quenches to Tf = Q we will argue in Sec. Ullthat the basic mechanism for the growth of L{t) can be properly seen as a 
progressive elimination of small domains with a faceted geometry, and that the zero temperature constraint, namely 
the unrealizability of activated moves, plays a crucial role. Elaborating on this we develop a scaling argument which 
allows us to determine analytically the behavior of several quantities. The results of this approach are compared 
in Sec. IIIII with the outcome of numerical simulations, providing a general agreement. In particular, for the total 
interface density p{t), which is related to the domains size by p{t)~^ <x. L{t), we find a power law behavior with z = 2 
in every dimension, as at finite temperature. The different character of the dynamics at Tf = Q is enlightened in 
Sec lII CI by considering the densities Pn{t) oc L{t)^^" of spins ai with a given degree of alignment n, this quantity 
being the difference between the number of aligned and that of anti-aligned neighbors of ai . These quantities provide 
information on the geometry of the interface and are shown to behave differently for quenches to finite Tf or to 
Tf = 0. While for shallow quenches, when the curvature driven mechanism is at work, one has /3„ = 1 for every n, 
for quenches to Ty = one finds n-dependent (and c?-dependent) values of /3„. For deep quenches with Tf > 0, a 
crossover is numerically observed fSec lIII B|) between an early stage (that can be rather long for small Tf) where the 
same behavior of quenches to Tf = is observed, to the late regime dominated by the usual curvature mechanism. 

Finally, we discuss the effects of temperature fiuctuations on the characteristic time tsc of the onset of scaling. Our 
numerical simulations (Sec. IIIip show that the behavior of tsc is very different in d = 2 and in d = 3. In d = 2, tsc 
is relatively small in quenches to Tj = and grows monotonously raising Tf. For Tf > and t < tsc one observes 
an approximate power law behavior with L{t) oc with an effective exponent Zcff{t) > 2, slowly converging 

to the asymptotic value. This explains why values of 1/z ~ 0.47 — 0.48 are often reported in the literature We 
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interpret the increase of tsc as due to the presence of the roughening length U{t, Tf) cx t^^^ competing with L(t) in the 
early regime. This interpretation is shown to agree with the results of numerical simulations. Moreover, we show how 
the effect of roughness can also be detected numerically in the behavior of h{x) (Eqs. (|3l4p ). Actually, over distances 
r < U(t,Tf) interfaces are not sharp, so that the Porod law ^ is not obeyed for x < Xfj{t) ~ U{t,Tf)/ L{t) ^ t~^/'^. 
For d = 3, instead, we find the opposite situation, tgc is very large at Tf — 0, while it is small for shallow quenches. In 
d~?>, U{t,Tf) grows at most logarithmically and hence is dominated by L{t) very soon causing no delays to scaling. 
Therefore, the mechanism leading to the increase of tsc when raising in d = 2 is not present here and, in shallow 
quenches, tsc is relatively small. Instead, when Tf = numerical simulations show a very long lasting transient. This 
is probably due to the constrained character of the kinetics where activated moves are forbidden. The very large value 
of tsc explains the anomalous values of 1/z ~ 0.33 — 0.37 sometimes reported in the literature [ll| for quenches to 
Tf 0. However, our simulations show unambiguously that l/zcff(t) is a growing function of i and, although at the 
longest simulated times it is still l/zeff{t) ~ 0.43, its behavior is consistent with an asymptotic value z — 2. 

This paper is organized as follows: In Sec|lT]we define the Ising model and develop scaling arguments to determine 
the behavior of p{t) and Pn{t) in quenches to Tf = 0. In Sec. HII] we discuss the data from simulations of quenches to 
Tf = (Sec. IIII A|) and to T/ > (Sec. IIIIB|) . showing the agreement with the results obtained in Sec. HH Sec. |IV] 
contains the conclusions. 



II. SCALING ARGUMENTS FOR QUENCHES TO Tf = 0. 

In the following we will consider the nearest neighbor Ising model described by the Hamiltonian 

H{[a]) = -JJ2<^^<^, (5) 

where cr^ = ±1 are the spin variables, are two nearest sites on a d-dimensional lattice and [a] is the configuration 
of all the spins. A purely relaxational dynamics without conservation of the order parameter can be defined by 
introducing single spin fiip transition rates w{ai —<Ti) obeying detailed balance. These quantities depend on Tf 
and on the local energy Ei — —Jui, rii being the degree of alignment, namely the difference between the number of 
the neighboring spins aligned with ai and that of the anti-aligned ones. Letting J — 1, transition rates are functions 
of Hi and Tf, namely w{ai — > —ai) = W{ni[a], Tf). For Tf = one has M^(rii[cr], Tf) = whenever > 0. 

In the remaining of this Section we will develop a scaling argument to determine the behavior of several quantities, 
among which L{t), in quenches to Tf = 0. 

We assume that the growing structure can be thought of as made of features, with a faceted geometry. Features 
are distortions of fiat interfaces or bubbles of spins. For the square lattice considered in the following, these are 
schematically drawn in Fig. [T] (upper part) in d = 2. 

A. Relation between the relaxation of a feature and the exponent z 

In order to have coarsening, features must be progressively removed [12i] . by fiipping all their spins. Let us define ti 
as the typical time to complete this process for a feature of size I. Our strategy is to relate n to the growth exponent 
z. In order to do that, let us notice that when, after a time t;, features of size I are removed, the typical scale of the 
system is increased by a quantity AI (x I, as shown in Fig. [1] (lower part). Assuming scaling, namely the presence of a 
single relevant lengthscale, the typical size I of a feature at time t must be of order L{t). Therefore dL{t)/dt ~ AI/ti. 
Let us anticipate what will be shown in the next Section, namely that t; oc P, with a = 2. Therefore we have 
dL{t)/dt (X L{t)^-°', and so L{t) cx t^/^ with 

z = a = 2. (6) 

In the following we consider the behavior of r; . 

B. Relaxation time of a feature 

We use the terminology of the case d = 2, for simplicity, but the argument is general. Let us consider the relaxation 
of an initially (at time t = 0) squared bubble, represented in Fig. [2l At zero temperature only spins with n < can 
be flipped. Therefore, referring to the situation of the upper part of Fig. ^ the first move is necessarily the flip of 
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FIG. 1: (Color online). Upper part: A flat surface in d = 2 with a distortion (left) and a bubble of reversed spins (right). 
Lower part: The increase of L{t) in d = 2 when features such as distortions (left) or bubbles (right) are removed. 



one of the four spins in the corners of the square. These moves trigger a sequence of successive flips, producing the 
shrinking of the bubble. Let us suppose, in order to simphfy the argument, that spins are flipped starting from the 
bottom of the box (actually the flipping of the spins proceeds on the average from each side, but this does not change 
our conclusions). Let us denote with hi (i = 1,2, ... ,1) the height of the i-th column of the box at time t, as shown 
in Fig. O Due to the zero temperature constraints, while the first and last column, with i = 1 and i = I, are always 
allowed to grow, due to the presence of the wall, all the other columns can do it only if at least one the nearest columns 




FIG. 2: (Color online). Relaxation of a bubble in d - 
typical configuration at a generic time. 



2. Upper: only the spins in the corner can be flipped initially. Lower: a 



is higher. Moreover, a column cannot decrease its height if it is not higher than at least one of the neighborings. 
With these rules, columns evolve until, at i = t;, all the spins in the box are flipped, the bubble disappears, and the 
process ends. 

Since with this dynamics an exact evaluation of t; is not possible, in the following we consider a slightly modified 
kinetics for which a determination of t; is allowed; we will then argue, checking this hypothesis numerically, that the 
modification of the dynamics does not change significantly the behavior of ti and, in particular, the exponent a. More 
precisely, we modify the original dynamics by introducing an additional constraint, namely — hi\ < 1. With this 
modification the problem can be mapped onto a diffusion equation for the variables hi. This result, which applies to 
the case d = 3 as well, is shown in Appendix I. For an interface described by a diffusion equation one has r; cx P, 
with a = 2. We argue that the same result applies to the original dynamics as well. The reason is the following: due 
to all the constraints discussed above, the heights hi are not independent, the typical differences \hi+i — hi\ do not 
grow very large and, for large I, they are independent of I. This is confirmed by looking at a simulation of the bubble 
shrinking. For large I, since the differences — hi\ are small as compared to the relevant scale I, we expect that 
the effect of the additional constraint does not change the exponent a. This last statement is convincingly confirmed 
by the results of numerical simulations, shown in Fig. [31 This figure refers to the simulation of a squared (cubic in 
d = 3) bubble, namely an Ising model on an P (l^ in d = 3) squared lattice with (say) up spins on the boundary and 
an initial condition of down spins in the interior. Averaging over several (10^ — 10^, depending on /) realizations of 
the thermal history, for each value of / we have computed ti as the time needed to revert the last spin. Fig. [3] shows 
that Ti oc with a = 2 is found with very good accuracy, regardless of dimensionality. This result confirms the 
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validity of the hypothesis according to which the exponent a is the same (a = 2) in the original Ising dynamics and 
in the modified kinetics considered in Appendix I. 

With this result Eq. ([6]) follows, namely z — 2 in every dimension. The exponent z is therefore the same as in 
quenches to finite temperatures. Assuming scaling, the size of domains L{t) is related by L(t) cx p{t)^^ to the total 
density of interfaces present in the system. The exponent z therefore gives informations on the number of interfaces, 
not on their geometry. In order to appreciate geometrical properties we will consider in the following other observables. 



C. Densities of spins with a given degree of alignment. 

Restricting again to the d = 2 case for simplicity ( the extension to the case d = 3 is straightforward and will 
be discussed in Appendix II). we introduce the density pn(/) of spins with a certain degree of alignment nt = n in 
a feature. In the following, we will only refer to interfacial spins: bulk spins with n = 4, whose behavior is trivial, 
will never be considered. Initially, in the squared bubble interfacial spins are those on the flat boundaries or in the 
corners, with n = 2 and n = respectively. During the evolution, as shown in the lower part of Fig. [2l all the possible 
values of n can be generated. The set of Pn{l) provides a geometric characterization of the interface. A representation 
of the typical geometry where a spin with a degree of alignment n occurs when the bubble shrinks is drawn in Fig. [4] 

While columns are growing, a generic profile of hi is made of steps, namely spins with n = 0, and flat parts with 
n = 2, as shown in Fig.[2l A step on site i can be randomly replaced by a flat part and the reverse is possible as well. 
As a consequence for sufficiently large values of I a finite fraction, independent of ^, of spins with n = and n — 2 
will be typically present. Simulations clearly show that these numbers, on average, do not depend on time (excluding, 
possibly, the initial and final stages of the process). According to this, the number of spins with n = or n = 2 is 
constant and proportional to the length of the interface, namely to I. Normalizing with the total number of spins 
we obtain P2(0 and po{l) oc l^^. The situation is very different for spins with n = —2 and n = —4. The 

former can only be produced when a last spin must be reversed in order to complete a row, an event happening, on 
average, every I moves. When this occurs, a single spin (out of P) with n = —2 is generated. Looking at a generic 
time, therefore, the typical density of such spins is p~2(0 ^ Finally, spins with n = —4 are only obtained when 
the last spin of the box must be reversed. In a bubble of P spins only one can be the last, and this happens once 
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FIG. 4: (Color online). Classification of interfacial spins in d = 2. The central spin is classified according to the degree of 
alignment. The dashed line is the typical shape of the interface associated to each type of spin during the bubble shrinkage. 



every P moves. Hence P-a{1) oc / ^. Again, scaling implies that we can identify I and L{t), leading to 

Pr,[L{t)]^L{t)-^-, (7) 

with 

/?2 = /3o = 1, (8) 

/3-2 = 3, (9) 

13-A = 4. (10) 
This argument can be extended to the case d — i (see Appendix II). The results are 

(ii = (i2=Po = 1, (11) 
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(i-2 = 3, (12) 

/3_4 = /3-6 = 4. (13) 

In the next Section we will compare these predictions with the outcome of numerical simulations, both in d = 2 and 
d = 3. 

III. NUMERICAL SIMULATIONS 
A. Quenches to T/ = 0. 

We have simulated systems of 2000^ and 576'^ spins in d = 2 and d = 3, respectively, on square (cubic) lattices. With 
these sizes, we have checked that finite size effects are not present in the range of times presented in the figures. The 
critical temperature of the model is Tc — 2.269 and Tc — 4.512 in d = 2 and d — i respectively (we set J — 1). Time 
is measured in montecarlo steps (mcs). In Fig. [5] we show the results for d — 2. Here we observe that a scaling regime, 
attested by the power law behavior of all the plotted quantities, sets in after a very short time tgc — 4 mcs. Best 
power law fits to the data (for t > 10) give P2 = 0.99 ± 0.02, /3o = 1.03 ± 0.03, /3_2 = 2.99 ± 0.04, /3_4 = 4.05 ± 0.05, 
and z = 1.99 ± 0.02. All the exponents /3„ and 2 are in excellent agreement with the determination made in the 
previous Section. 

When scaling holds, according to Eqs. (HI,® the equal time correlation function, behaves as 

Gagir,t,t)^hir/L{t)), (14) 

with h obeying Porod law ^ in the case of sharp interfaces. In the case considered here, Gag{r, t, t) = G{r, t, t), since 
Gst{r,t,t) = a.t Tf — 0. According to Eq. p^ . when scaling holds curves of Gag{r,t,t) for different times should 
collapse when plotted against r/L{t). This is observed in Fig. [H] (left part). Besides, the Porod law Q is very neatly 
obeyed. 

In Fig. © a plot analogous to that of Fig. O is made for d = 3. After a time around 10 mcs a power-law behavior 
sets in for all the pn[L{t)\. Fitting the curves with power laws we find a residual time dependence of the exponents /3„, 
since their value changes measuring them in different timewindows. This is particularly evident for /3_2- This indicates 
that preasymptotic corrections to scaling are not completely negligeable in the timedomain of our simulations. Best 
power-law fits for t > 10 yield (3^ = 0.97 ± 0.04, (32 = 1.11 ± 0.06, /3o = 1.03 ± 0.04, [3-2 = 3.3 ± 0.1, /3_4 = 3.8 ± 0.1, 
and (3-Q = 3.9 ± 0.1. Taking into account the presence of preasymptotic corrections we regard these values as being 
consistent with the results of the previous Section. Regarding L{t), instead, the data (in the inset) do not show 
a satisfactory power law. The curve is bending upwards on the double logarithmic plot and the exponent is not 
consistent with the expected value z = 2. In order to clarify this point we have computed the effective exponent 

1 dlnLit) , , 

^ ^ (15) 



Zeff{t) dint ' 

which is shown in the inset of Fig. [SI In an early stage, when scaling does not hold, this quantity grows exponentially 
as described by linear theories 13]. Then, after reaching a minimum of order 0.3, l/ze//(t) it keeps slowly, but 
steadily, increasing. Its value measured at the longest times is around l/zeff — 0.43. Actually, to the best of our 
knowledge, the expected exponent z = 2 has never been reported. Previous simulations on much shorter timescales 
observed [ll| an exponent of order 1/3 and sometimes the very existence of dynamical scaling has been questioned. 
Notice that the value 1/3 is comparable to the value 1/zeff — 0.3 of the effective exponent around its minimum, 
a fact which may explain what reported in [ll[. Our data are consistent with a possible asymptotic value z = 2 
although much larger simulation efforts would be needed for a definitive evidence. In any case, the data show that 
preasymptotic effects are quite relevant, and exclude that a well defined exponent can be measured, up to i ~ 2 • 10^ 
mcs (for longer times, not reported in the Figure, finite size effects are observed). A rough extrapolation suggests 
that times at least 10 times larger (t ~ 2 x 10^ mcs) are needed to observe l/zeff{t) ~ 1/z ~ 1/2 (meaning a lattice 
size of order (2 • 10^)^ in order to be finite size effects free). 



B. Quenches to Tf > 0. 



When the quench is made to a finite final temperature all the constraints imposed by T/ = are removed. In the 
quench to Tf = 0, as seen in Sec. IIII Al spins with n > are present in the system, but they cannot be updated 
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FIG. 5: (Color online). Quench to T/ = in d = 2. The densities pn[L{t)] are plotted against L{t) — p{t) ^ . In the inset L{t) 
is plotted against time. 



because this would increase the energy. When Tf > also these can be updated, although with a small probability 
for small Tf. For shallow quenches (T < Tc) the typical times t„(T/) = W{n.i[a],Tf)^^ associated to microscopic 
moves are small, and, in particular, much smaller then the timescales over which the non-equilibrium behavior of 
interfaces takes place. Therefore we expect that during phase-ordering an interface be in quasi equilibrium, namely 
it will have the same values of pn[L{t)] of an equilibrium interface [14,] of length I = L{t) at the same temperature 
of the quench. In order to check this conjecture we have performed the following simulations: In c? = 2 we have 
prepared an Ising system with a spanning vertical interface in the middle and antiperiodic boundary conditions in the 
horizontal direction; subsequently, we let it evolve at a constant temperature T and, during the evolution, we have 
measured the densities Pn(0- These quantities are compared to the quantities pn[L{t)] measured in a quench to the 
final temperature Tj — T. In both the situations we have implemented a fast dynamics where spins with n = 4 are 
not allowed to flip. This no bulk flip (NBF) dynamics has been frequently used in the literature [H, Apart from 
its numerical efficiency, it has the advantage of isolating the aging behavior of the system. The reason is that since, 
as already pointed out in Sec. HI the stationary terms in Eq. ([T]) are produced by the flipping of the spins inside the 
domains, by preventing bulk moves one is left only with the dynamics of the interfaces which is responsible for the 
aging term of Eq. ([1]). On the other hand, it has shown [IG'I also that the NBF rule does not change the properties 
of the aging terms in the large time domain. In the simulation of the single interface we use this no bulk flip (NBF) 
dynamics also as a tool to maintain a single interface in the system at all times. With the standard dynamics spins 
can be reversed in the bulk, creating additional interfaces that may interact with the original spanning interface. 
On the other hand, with NBF dynamics, the spanning interface remains unique and well defined at all times. In 
order to avoid the complications arising when comparing Pn{l) with pn[L{t)\ in the two kind of simulations, which 
would require the comparison of the size I of the interface in the equilibrium simulation with the length L(t) in the 
corresponding quenched system, in Fig[7|we have plotted the ratios between different /?„. Since these quantities do 
not depend on I or on L{t) respectively in the two cases, they can be directly compared. After a brief transient, the 
single interface (main figure) reaches the stationary state and the ratios between different pn take time-independent 
values. The same is true for the quenched system, shown in the inset. We find that the asymptotic value of the ratios 
is the same with good accuracy in the two systems. This confirms our claim that in a shallow quench the values of 
the densities pn[L{t)] are equal to the corresponding quantities Pn(l) in an equilibrium interface of size Since the 
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FIG. 6: (Color online). Quench to Tf — in d = 3. The densities p„[L{t)] are plotted against L{t) — p{t) ^. In the insets L{t) 
and l/zeff{t) are plotted against time. 



latter are finite and time-independent this implies that in a shallow quench /7„[L(t)] cx p{t). Therefore, instead of 
Eqs. P fTO]) . and Eqs. (|imi3p one must have 

Pn = 1, (16) 

for every value of n. This is shown to be true in Fig. [8] for systems in d = 2, 3 quenched to Tf — 2. Best power-law 
fits (for t > 10) yield (32 = 1.00 ± 0.01, (3o = 1.00 ± 0.01, (3-2 = 1-00 ± 0.01, /3_4 = 1-01 ± 0.1, in d = 2 and 
Pi = 1-00 ± 0.01,/32 = 1.00 ± 0.01,/3o = 1.00 ± 0.01, /3_2 = 1.00 ± 0.01,/3_4 = 1.01 ± 0.01, P2 = 1.01 ± 0.02, in d = 3. 
Notice that the values of /3„ with n < are very different from the case with Tf — and can be used to distinguish 
the two kinds of dynamics. Regarding the value of the exponent z, the curvature driven mechanism implies z — 2. 
Actually, this is fomid with very good accuracy in d = 3 (we find l/z — 0.502 ± 0.04 in the range t G [10^ — 10^]). 
Instead, for d = 2 one observes a slightly larger exponent, since 1/z is of order 0.48 in the region of the largest 
simulated times. We will comment later on this point. 

Since the dynamics is diff'erent in shallow quenches or in quenches to = we expect to see a crossover phenomenon 
at intermediate temperatures. Namely, for every > a crossover time should exist separating an early stage where 
the dynamics is of the Tf = type, with /3„ given in Eqs. ([8l fT0|) for d = 2 or in Eqs. (jlHllSp for d = 3, from a late 
stage where the finite temperature scalings ([T5)) set in. For a class of spins with a given degree of alignment n the 
crossover between the early and the late kind of dynamics occurs when spins with the considered n start to be created 
by means of activated moves. The crossover time, therefore, should be of order t„(T/) ~ W{n, Tf)~^, and is therefore 
different for spins with different n. At the crossover time t„(T/) a typical crossover length 

L„(T/) ^ r„(T/)i/^ = W{n,Tf)-^/^ (17) 

is associated. In Figs. 191101 the pattern of crossover described above can be observed. Here we see that, practically for 
all the Tf considered, the behavior of the densities Pn{t) is initially analogous to that of the Tf = case (Eqs. ((8] fT0)) 
or Eqs. (|lim3p in d = 2 or d = 3 respectively). This regime lasts until L{t) ~ Ln{Tf), where Pn{t) start to behave as 
in shallow quenches fEg. [TB|) . For very small Tf, Ln{Tf) is outside the range of simulated times (this explains why, 
for instance, the curves with T/ = and Tf = 0.25 can be hardly distinguished in d = 2). Increasing Tf gradually, 
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FIG. 7: (Color online). The ratios between p_2(/), po(l), and p2(l) are plotted against t for an equilibrium interface in d = 2 
a.t T — 2 (main) and in a quench to T/ = 2 (inset). NBF dynamics is used in both cases. 



Ln{Tf), whose values obtained from Eq. (|17p are marked with vertical segments across the curves (when within the 
simulated times), become progressively smaller. One observes that the crossover phenomenon occurs at different times 
for spins with different n, and that the estimate ^T7\ agrees reasonably well with what observed. 

Let us now come back to the value of the exponent 1/z in d = 2, which, as already observed regarding Fig. [51 has a 
value slightly smaller than the expected one 1/z = 0.5. In order to make more precise statements we have measured 
the effective exponent, which is shown in Fig. [TT]for various temperatures. For T — the effective exponent initially 
rise to a maximum for the reasons already discussed for the case with d = 3. Then it goes down to a minimum and, 
later, reaches the asymptotic value 0.5 already at times of order t ~ 300 mcs. As Tf is increased the pattern is similar 
but the initial minimum is depressed and delayed so that for the largest temperatures considered l/zeff(t) has not yet 
reached the asymptotic value at the longest simulated times. Although the expected final value l/ze// — ^/z = 1/2 
is not in doubt, a rough determination of this exponent in a simulation may lead to a smaller value, as sometimes 
reported |17!]. This behavior can be interpreted as due to the presence, beside L{t), of another length, the roughness 
of the interfaces. Equilibrium interfaces are rough for T > T^. The roughening temperature Tji vanishes for d = 2 
while < Tfi < for c? = 3. An interface spanning a box of linear size I in equilibrium at the temperature T has a 
typical width ui{T) given by [l^ 

f a2(T)yfl for d = 2 
uiiT) ^ I a:i{T) for d = 3,T<TR (18) 
[a3(r)ln/ for d = 3,T>TR. 

In the phase-ordering kinetics it has been conjectured by Villain [l^ that the role of I in Eq. (fT8|) is played by L{t). 
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FIG. 9: (Color online), d = 2. p-2[L{t)] (left) and p-4,[L{t)] (right) are plotted against L{t) = p(t)~^ (NBF). Vertical segments 
across the curves represent the crossover lengths L-2{Tf) and L-4,{Tf) of Eq. (when reached in the simulation). 
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FIG. 10: (Color online), d = 3. p_2[L(t)], p-i[L[t)] and p-6\L{t)] are plotted against L{t) = pit)''^ (NBF). Vertical segments 
across the curves represent the crossover lengths L_2(T/), L-i{Tf) and L-^{Tf) of Eq. (|17|) (when reached in the simulation). 



The non-equilibrium width U{t, Tf) should than behave as 

r a2{Tf),/L{t) 
U{t,Tf) cx <^ asiTf) for 
[ a,iTf)lnLit) 



for d^2 

d = 3, Tf < Tr (19) 
for d^3,Tf > Tr. 



According to these expressions, in the large time limit U{t, Tf) can always be neglected with respect to L{t). However, 
there can be an initial regime, for t < tsc, where U{t,Tf) produces a correction to scaling. In this range of times 
we expect a (time dependent) effective exponent Zeff(t) ^ 2 to be observed. Given the behaviors p^ . tsc may be 
sufficiently large to produce observable effects for d = 2 while we expect it to be too small to significantly affect the 
scaling behavior for d = 3. Actually, we have already observed (see Fig. [8]) that, differently from d = 2, in d = 3 the 
effective exponent quickly converges to the value z = 2 for quenches to Q < Tf < T^. According to our hypothesis, 
since a(Tf) is an increasing function of Tf, while L{t) is roughly temperature independent, the convergence towards 
the asymptotic z = 2 should be delayed increasing Tf: This is actually observed in Fig. 1111 In order to check further 
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FIG. 11: (Color online). The effective exponent l/zeff{t) for quenches in d = 2 (NBF). 



the consistency of this hypothesis we have computed the temperature dependence of a(Z, Tf). From a set of simulations 
of a single interface as those described above in this section we have extracted the equilibrium width of the interface 

as ui{T) = {^jY^'j^i{^/l)[xj{t) — 1/2]'^), where j is the vertical coordinate in the simulation box, Xj{t) is horizontal 

coordinate of the interface position at a generic time t, and (• • •) is an average over thermal realizations. Since in 
equilibrium ui{T) does not depend on time we have also averaged ui{T) over time in order to reduce the noise. Fig. 1121 
shows that the behavior is obeyed for I sufficiently large (the larger the lower is Tj). Extracting 02 (Ty) we find 
a linear relation 

a2{Tf)^ATf, (20) 

where A is a constant {A ~ 0.32). We can evaluate tsc from the condition L{tsc) = U{tsc,Tf). In a quench from 
high temperature L{t) start growing from an initial value L{0) ~ 1. For low temperatures, since 02(1/) is very small, 
L(0) is larger than U{0,Tf) and hence U{t,Tf) is negligible from the beginning. In this case scaling can set in very 
early, after the microscopic time i* ~ 1 necessary for the formation of domains of the equilibrium phases which is 
practically independent of Tj. For larger temperatures there is a transient during which U {t, Tf) cannot be neglected. 
Using Eqs. (|19l20p one obtains L{tsc) = A^Tj. Since Lit) ^ t^/'^ is roughly obeyed also for t < tsc (the effective 
exponent is always in the range [0.45 — 0.5]) one can estimate tsc « T"'. In Fig. [13] we have plotted L(tsc) for different 
values of Tf. This quantity have been obtained as follows: From the data of Fig. [TT]we have estimated tsc as the 
time when l/zeff{t) reaches the value 0.48 (clearly, we refer to the asymptotic increase of l/zeff{t), for t > 100, 
not to the early maximum). Successively, from the numerical data for L(t) we have extracted L(tsc)- The picture 
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shows agreement with the prediction of our hypothesis, namely a constant vahie of L{tsc) at low temperatures and a 
behavior L(isc) ~ for larger temperatures. 

The interplay between U {t, Tf) and L{t) can also be observed in the behavior of the equal time correlation function, 
which, when scaling holds, should behave as in Eq. (|14p . with h obeying Porod law ([4]) in the case of sharp interfaces. 
However, as already discussed, the presence of U{t,Tf) introduces a correction to scaling in an early regime when 
L(t) has not yet grown sufhciently larger than U(t,Tf). Moreover, due to roughness, interfaces are not sharp. Then, 
both scaling and the Porod law are expected to be violated for r < U{t,Tf), that is for x — r/L(t) ^ Xfi{t) — 
U{t,Tf)/ L{t) — a2{Tf)L^^/^{t), namely in a range of x that shrinks in time but that may be appreciable for large 
Tf. Actually this is observed in Fig. [TH While curves of Gag{r,t,t) for different times collapse when plotted against 
r/L(t) for the quench to Tf = 0, as discussed in Sec. IIII Al and the Porod law is also verified, when the quench to 
Tf > is considered one observes significant scaling violations in the region of small x. In this regime, the curve 
definitely deviates from the linear Porod law. As time goes on, these violations become weaker and the curves seem 
to approach the same behavior as for Tf = 0. For intermediate temperatures similar, but less pronounced, violations 
are also observed. 
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FIG. 13: (Color online). L{tsc) is plotted against Tf 



IV. SUMMARY AND CONCLUSIONS 



In this paper we have investigated the role of Tf in the phase-ordering kinetics of the Ising model with single spin 
flip dynamics. At T/ = the dynamics is characterized by faceted interfaces and by the kinetic constraint of the 
impossibility of activated moves. At > interfaces are curved and rough (for T > Tn in d = 3). We have shown 
that, while the exponent z regulating the decay of the total density of interfacial spins is not changed by the different 
geometry of interfaces in quenches to Tf — or Tf > 0, other quantities, such as the exponents /?„ describing the 
behavior of the densities of particular classes of spins do change. The existence of two different dynamical mechanisms 
induces a crossover pattern for finite Tf. Besides, in d = 2 for T/ > the roughening length competes with L(t) in an 
early stage, delaying the realization of dynamical scaling, as it is evidenced by the time dependence of the effective 
exponent z^ff and by the breakdown of the Porod law at small r/L{t). 

This whole pattern of behaviors is due to equilibrium properties of the interfaces. Therefore we expect to observe 
similar behaviors for the Ising model with conserved dynamics. 



APPENDIX I 
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FIG. 14: (Color online). G{f,t,t) is plotted against x — r/L{t) for a quench to T = (left) or T = 2 (right). Different curves 
correspond to several times (the same in the two figures) between t = IQ and t = 10* (from bottom to top, for T = 2). The 
dot-dashed line is the Porod law y = 1 — ax, where a (the same in the two pictures) is obtained as the best fit of G{f,t,t) at 
T = 0. Points are joined by a piecewise continuous line as a guide for the eye (for clarity the first two points, the one in the 
origin and the following, are not joined). 



We consider a generic profile of the hi and denote with {h} this configuration. Let us introduce the probability 
P{{h}, t) of having such a configuration at time t, and the conditional probability P{{h'}, t'\{h}, t) of having {h'} at 
time t' provided that the configuration {h} was found a,t t < t' . One has 



{h,{t + At) ~ h,{t)) ^ {h'^~h,)P{{h'},t + At\{h},t)P{{h},t). 

{h}.{h'} 

From the master equation, the conditional probability can be written as [20| 



(21) 



pah'},t+At\{h},t) = 



3 h'^'^hj 



d[h'}.{h} + ^tJ2 E w{h'^ ^h'^)S[h"}.{h}+0{At^) (22) 



where we have introduced the transition rate w{hj — > h'j) for moving the height hj of the j-th column to h'j and 
S{h"},{h} is the Kronecker function between configurations {h"} and {h}. We have assumed that the w are single spin 
flip transition rates and hence transitions between configurations can be obtained by summing over j. Apart from 
this, the w are still generic (not necessary those of the zero temperature Ising model) at this stage. Inserting Eq. 
into Eq. (j2ip . due to the term (h^ — hi) only contributions with {h'} ^ {h} and j = i survive, and one has 



(/i,(t + Ai)-/i,(t)) = Ai J2 {K-h,)P{{h},t)w{h,^h'i). 

{h},{h']^{h) 

Since in an elementary move /i^ — > /i^ ± 1 is only allowed, introducing m = ±1 the last Equation reads 

{K{t + At) - h,{t)) ^ AtY^Y^ 

mw(hi ^ hi + m)P{{h}, t). 

{h} m 



(23) 



(24) 
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Taking the continuum limit At — > yields 

^M)l^J2Y,rawih,^h, + m)Pi{h},t). (25) 

{h} rn 

This equation has been obtained without any approximation. We now have to specify the form of the w{hi hi + m) 
in order to reproduces the original rules of the zero temperature Ising model. However, in order to have a tractable 
model, we consider transition rates which correspond to the Ising model with the additional constraint 

\h,+i{t) - h,{t)\ < 1 Vz, (26) 
as discussed in Sec. Ill Bl Starting with the case d = 2, we define them as 

wih, -^h, + m)^ ^V^/i, + F.t\{h}), (27) 
where V^/i; — /i^-i + /li+i — 2hi is the discrete Laplacian in one dimension, and 

In order to see this let us notice first that if the configuration {h} satisfies the condition (pS)) . with the transition 
rates (|27p that constraint will never be violated by the later evolution. In fact, let us focus on site i and suppose 
that the transition hi —^ hi + 1 is going to be attempted. This transition would violate Eq. (j26p if /i^+i = hi — I. In 
this case, however, it easy to check that for every consistent with the condition (|26p it is ■w{hi hi + 1) = 0. 
The same argument can be repeated for every configuration. Having proved that (|26|) is fulfilled by the transition 
rates ^17} we can restrict ourselves to consider the only cases allowed, namely those with V^/i^ = 0,±1,±2. When 
V^/ii = the only configuration where the move hi ^ hiztl could be attempted without violating the constraint (|26|) 
is that with hi^i ^ hi ^ /li+i. In this case, in the original Ising model the move is forbidden, which agrees with 
Eq. ([77]) giving w{hi hi + m) = in this case. Coming to the cases with V^hi = ±1 (which are realized, for 
instance, when hi^i — hi — hi^i ^ 1), in the original Ising model moves with m sign (V^/i^) < are energetically 
forbidden, while those with m sign {V^hi) > occur with a probability 1/2, using Glauber transition rates. This 
agrees with Eq. ([27|l . Analogously, when V^/ii = ±2 (when, for instance, hi+i = hi^i and hi = hi+i T l)j in the 
original Ising model moves with m sign (y'^hi^ < are forbidden, while those with m sign (V^/ii) > lower the 
energy and occur with a probability 1, providing again agreement with Eq. (|27p . 

Let us now insert the transition rates ([27|) into the evolution Equation (f25|) . obtaining 



d{h,{t)) 



E E ^^'h^pm, + E E m)pm, t). (29) 



dt 

{h} ™ {h} 



Performing the sum over m one has 

d{h,{t)) 



= (V^/i.) + {Ft) (30) 



dt 

where Ft{{h}) = Fj;^\{h}) — F^~^\{h}) . Let us now consider the last term on the r.h.s. of Eq ([50]) . We want to show 
that it can be neglected. We will show it separately, for all the possible values of V^hi, namely V^/i^ = 0, ±1, ±2. For 
V^hi — 0, from Eq. (|28p . it is Ft[{h}) = 0. Let us consider now the contributions with V^ft.^ = ±1. This situation 
corresponds to spins with n = 0, or steps in the terminology of Sec. Ill CI These that can be flipped from = 1 to 
tJ = — 1 and back without energy costs. Therefore the two values of the spin in this case occur with equal probability 
(i.e. 1/2) and, as V^hi changes its sign when the spin is reversed the contributions V^hi = ±1 cancel in (F^)- The 
case V^/ii = —2 can never be realized in the kinetics because it would require a move with energy increase. Therefore, 
the terms with V^hi = 0, ±1, —2 do not contribute to the r.h.s. of Eq. ((30|) . This is no longer true for V^/ii = 2. This 
term corresponds, in the language of Sec. Ill CI to a spin with n = —2 which, as explained in Sec. Ill CI are created 
only when the last spin has to be reversed to complete a row. On average this happens once every I moves. Therefore, 
although the contributions with V^/i^ = 2 do not strictly vanish, they provide a contribution (F^) oc l/l, which can 
be neglected in the large-Z limit. Then we arrive at the diffusion equation 



(31) 
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1 V^h, 1 










1 


-1/4 


2 


r -1/2 for Vl=OorVl = 




\ else 


3 


1/4 


4 






TABLE I: The values of the function fi({h}). i^y) is the discrete second derivative along x (y). 



Let us turn now to the case d = 3. Similarly to the case d = 2, it is easy to check that the transition rates 



with the following form of Fj;"^\{h}) 



h, + m)^jV'h, + Fl"'\{h}), 



(32) 



Pi"^ {{h}) — ^|V^/l,;|(5,„sjg„(v2;i,),_i + fii{h})6,nsign{V^h,), 



(33) 



where W^hi is the discretized laplacian in d = 3 and with the values of fi{{h}) given in table HI satisfy the con- 
straint (|26|) and reproduce the Glauber transition rates of the original Ising model at T = 0. Proceeding as for d = 2 
one arrives at 



(34) 



By reasoning as in d = 2, one concludes that only sites with |V^ft.j;| — 1,2,3,4, corresponding to spins with Ei ^ 
contribute to (F^), but these can be neglected for large I. Hence one arrives also in this case to a diffusion equation 



d{h,{t)) 
dt 



(35) 



APPENDIX II 



Let us consider the shrinkage of a cubic bubble of linear size I, and extend the argument developed in Sec. Ill CI to 
the case d = 3. The interfacial spins can be classified according to n, as shown in Fig. 1151 Suppose again that the 
interface grows from the bottom of the bubble. Let us define hi as the height of the z— th column, with i — I, 
running on the two-dimensional lattice. While columns are growing, the profile of hi is made of fiat parts (spins with 
n = 4), edges (spins with n=2) and corners (spins with n = 0). We make again the hypothesis that the probability 
of their occurrence is finite and constant. By reasoning analogously to the d = 2 case, since these spins belong to 
the growing surface, their number is proportional to P. Hence Pi{l) oc p2{l) oc Pq{1) oc l^^ . Let us imagine to reverse 
the spins from the bottom, level by level. A spin with n — —4 is only produced when the last spin of a certain level 
has to be reversed. All the I levels are completed in a time r/ (x so that there is a number proportional to of 
spins with n = —4 in a unit time. This implies p^a{1) oc Z^*. Spins with n = —2 are generated analogously to the 
spins with n — —2 in d = 2, namely when the last spin must be reversed in order to complete a row. rows must be 
completed in a time r oc Z^ in order to reverse all the spins of the bubble. Therefore there are P spins with n = —2 in 
a unit time. In conclusion p-2{l) oc Finally, spins with n = —6 are only formed when a growing column reaches 
the top of the bubble, with the condition that all the nearest column have already reached the top. Due to the kinetic 
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FIG. 15: (Color online). Classification of interfacial spins. 



constraints, these neighboring column must themselves have at least one nearest column of equal or higher eight. One 
can iterate this argument until a border is reached. Since the border is a distance of order I away, one concludes 
that every I spins only one can have n = — 6. Therefore a number proportional to I of spins are generated when the 
columns reach the top, and this event happens once every r; cx P moves. Then p-&{l) cx Assuming scaling, and 
identifying / with L{t) in the phase-ordering kinetics one arrives at Eqs. (|llj|13p . 
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